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X chromosomes are preferentially transmitted tlirough females, wliicli may favor tlie accumulation of X-liniced alleles/ 
genes with female-beneficial effects. Numerous studies have shown that genes with sex-biased expression are under- or 
over-represented on the X chromosomes of a wide variety of organisms. The patterns, however, vary between different 
animal species, and the causes of these differences are unresolved. Additionally, genes with sex-biased expression tend to 
be narrowly expressed in a limited number of tissues, and narrowly expressed genes are also non-randomly X-linked in 
a taxon-specific manner. It is therefore unclear whether the unique gene content of the X chromosome is the result of 
selection on genes with sex-biased expression, narrowly expressed genes, or some combination of the two. To address this 
problem, we measured sex-biased expression in multiple Drosophila species and at different developmental time points. 
These data were combined with available expression measurements from Drosophila melanogaster and mouse to reconcile the 
inconsistencies in X-chromosome content among taxa. Our results suggest that most of the differences between Drosophila 
and mammals are confounded by disparate data collection /analysis approaches as well as the correlation between sex bias 
and expression breadth. Both the Drosophila and mouse X chromosomes harbor an excess of genes with female-biased 
expression after controlling for the confounding factors, suggesting that the asymmetrical transmission of the X chro- 
mosome favors the accumulation of female-beneficial mutations in X-linked genes. However, some taxon-specific patterns 
remain, and we provide evidence that these are in part a consequence of constraints imposed by the dosage compensation 
mechanism in Drosophila. 

[Supplemental material is available for this article.] 



Sexual dimorphism is widespread among animal species. Because 
males and females share nearly identical genomes, differences in 
phenotjfpes are in part the result of the biased expression of certain 
genes in one sex or the other (EUegren and Parsch 2007); these 
genes are said to have "sex-biased expression." Selection pressures 
can differ between males and females, and these sex-speclflc or 
sexually antagonistic selection pressures are expected to act on 
genes with sex-biased expression (Vicoso and Charlesworth 2006; 
Innocenti and Morrow 2010). X chromosomes are found in fe- 
males two-thirds of the time, whereas autosomes are equally 
transmitted via the two sexes. The net selection pressure acting on 
X-linked genetic Vciriation may therefore be biased in favor of fe- 
male-beneficial alleles (Rice 1984; Connallon and Clark 2010; Fry 
2010). Assuming that genes with sex-biased expression perform 
sex-specific-beneficial functions (Innocenti and Morrow 2010), 
the X-llnkage of genes with sex-biased expression can thus be 
a valuable readout of historical sex-specific selection pressures 
(Vicoso and Charlesworth 2006). 

The degree to which the sexes are dimorphic also varies across 
tissues. For example, some somatic tissues are highly homoge- 
neous between the sexes, other somatic tissues are unique to one 
sex (e.g., male accessory glands in insects, the prostate in mammals), 
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and the gonads are extremely differentiated between males 
and females. Each tissue expresses a subset of the genes in a 
genome — the more sexually dimorphic the tissue, the more dif- 
ferentiated is the pattern of female versus male gene expression in 
that tissue Oin et al. 2001; Arbeitman et al. 2002, 2004; Parisl et al. 
2003, 2004). In addition, the expression breadth of Individual 
genes can range from broad (i.e., the gene is expressed in multiple 
tissues/developmental stages/environments) to narrow (expressed 
in very few tissues/stages/environments). Sex-specific selection 
pressures may differ if a gene is narrowly expressed in a sex-limited 
reproductive tissue or a tissue shared by both sexes (Meisel 2011). 
Furthermore, sexually antagonistic selection pressures are expec- 
ted to be greatest in sexually mature adults and minimal in sexually 
immature juveniles (Chippindale et al. 2001; Rice and Chippindale 
2001; Gibson et al. 2002). The fitness of a genotype is a composite 
of the performance of that genotype integrated across all tissues, 
developmental stages, and environments in which the gene(s) is/ 
are expressed. Therefore, if one is to interpret the selection pres- 
sures acting on genes based on their X-Iinkage, one must consider 
both when and where each gene is expressed in the organism. 

Multiple experiments have revealed that genes with male- or 
testis-blased expression are under-represented on the Drosophila X 
chromosome, while there is an excess of X-linked genes with fe- 
male- cmd ovary-biased expression (e.g., Paris! et cil. 2003; Ranz 
et al. 2003; Sturgill et al. 2007; Mikhaylova and Nurminsky 2011). 
Additionally, the right arm (XR or neo-X) of the Drosophila pseu- 
doobscura X chromosome (which is autosomal in most other 
species) also contains a significant deficiency of genes with male- 
biased expression (Sturgill et al. 2007), suggesting that the paucity 
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of X-linked genes with male-biased expression is an emergent 
property of the Drosophila X chromosome cind not Intrinsic to 
the proto-X. 

Three hypotheses have been put forth to explain the de- 
ficiency of genes with male-biased expression on the Drosophila X 
chromosome. First, sexually antagonistic selection may prevent 
the accumulation of X-Iinked male-beneficial mutations because 
the X chromosome is preferentially transmitted through females 
(Rice 1984; Connallon and Clark 2010). Second, male meiotic sex- 
chromosome inactivation (MSCl) may prevent the expression of 
X-linked genes during spermatogenesis (Betran et al. 2002; Hense 
et al. 2007; Vibranovski et al. 2009), which could decrease the 
proportion of X-linked genes with testis-biased expression. Third, 
dosage compensation, via hyperexpression, of the hemizygous X 
chromosome in males may limit the degree to which X-linked 
genes can acquire male-biased expression (Vicoso and Charlesworth 
2009; Bachtrog et al. 2010). It is still unclear, however, the extent to 
which sexual antagonism, MSCI, and "dosage limits" each con- 
tributes to the deficiency of X-linked genes with male-biased 
expression. 

In mammals, the X chromosome also harbors a unique col- 
lection of gene types, but in a way that is thought to differ from the 
patterns in Drosophila (Vicoso and Charlesworth 2006; Gurbich 
and Bachtrog 2008). For example, the Drosophila accessory gland 
and the mammalian prostate are both male-limited reproductive 
organs that contribute seminal proteins to the ejaculate. The 
Drosophila melanogaster X chromosome is depauperate in genes 
encoding accessory gland proteins (ACPs) (Wolfner et al. 1997; 
Swanson et al. 2001; Mueller et al. 2005) and genes with accessory- 
gland-biased expression (Mikhaylova and Nurminsky 201 1), while 
there is an excess of X-linked human genes with prostate-biased 
expression (Lercher et al. 2003). 

In addition, unlike in Drosophila, the mammalian X chro- 
mosome has an excess of genes with testis-biased expression 
(Wang et al. 2001; Lercher et al. 2003; Mueller et al. 2008). A closer 
examination, however, has revealed that this excess is limited to 
genes with pre-meiotic expression, whereas post-meiotically 
expressed genes are under-represented on the X chromosome be- 
cause of MSCI (Khil et al. 2004). Similarly the paucity of X-Iinked 
testis-expressed D. melanogaster genes may be limited to genes 
expressed during later stages of spermatogenesis such as meiosis 
(Vibranovski et al. 2009; Gan et al. 2010). Furthermore, new genes 
with testis- or male-biased expression tend to be preferentially lo- 
cated on the X or Z chromosome in Drosophila (Zhang et al. 2010a), 
mammals (Zhang et al. 2010b), emd chicken (Ellegren 2011). 
Therefore, the X- (or Z-) linkage of testis-expressed genes may not 
differ substantially between taxa. 

Analyses of non-reproductive tissues have also revealed sim- 
ilarities in the X-linkage of genes with sex-biased expression in 
Drosophila and mammals. For example, genes with male-biased 
expression in mouse brain are over-represented on the X chro- 
mosome (Yang et al. 2006), as are genes with male-biased expres- 
sion in the D. melanogaster head (Chang et al. 2011). Also, genes 
with female-biased expression in brain, liver, and muscle are over- 
represented on the mouse X chromosome (Yang et al. 2006), and 
there is evidence for an excess of X-Iinked D. melanogaster genes 
with female- or ovary-biased expression (Parisi et cil. 2003; Ranz 
et al. 2003; Mikhaylova and Nurminsky 2011). These results sug- 
gest that, while there are some differences in gene content between 
the Drosophila and mammalian X chromosomes, they may be 
more similar than previously thought (Vicoso and Charlesworth 
2006; Gurbich and Bachtrog 2008). The remaining differences 



between the X chromosomes could be the result of differences in 
dosage compensation (Vicoso and Chcirlesworth 2006), different 
dominance coefficients of sexually antagonistic genetic variants 
(Rice 1984), or artifacts of how expression data are collected and 
analyzed in each taxon. To distinguish between these hypoth- 
eses, we must analyze comparable data from Drosophila and mouse 
using the same approaches. 

In addition to the X chromosome harboring a non-random 
proportion of genes with sex-biased expression, the expression 
breadth of X-linked genes also differs from that of autosomal 
genes: The D. melanogaster X chiomosome contains a deficiency of 
narrowly expressed genes (Mikhaylova and Nurminsky 2011), 
while there is an excess of X-linked narrowly expressed human 
genes (Lercher et al. 2003). Genes with sex-biased expression in 
non-reproductive tissues tend to be more narrowly expressed in 
both D. melanogaster and mouse (Mank et al. 2008; Meisel 2011). 
Additionally, D. melanogaster genes with female-biased expression 
are often highly expressed in ovary, and ovary-expressed genes 
tend to be broadly expressed (Parisi et al. 2004; Meisel 2011). It 
is therefore not clear whether the selection pressures that drive 
differences in gene content between the X chromosome and 
the autosomes are acting on genes with sex-biased expression, 
narrowly expressed genes, or some combination of the two 
(Mikhaylova and Nurminsky 2011). 

The X-linkage of genes with sex-biased expression is an in- 
formative metric of sex-specific selection pressures (Vicoso and 
Charlesworth 2006), but to properly interpret these patterns and 
compare across taxa, one must control for when and where genes 
are expressed. To address the limitations of previous studies, we 
measured sex-biased expression in the heads of four Drosophila 
species (D. melanogaster, D.pseudoobscura, Drosophilawillistoni, and 
Drosophila mojavensis) using a combination of microarrays and 
high-throughput RNA sequencing (RNA-seq), in whole flies from 
three species (D. melanogaster, D. pseudoobscura, and D. mojavensis) 
using RNA-seq, in thorax and abdomen of D. willistoni using RNA- 
seq, and in larvae of D. melanogaster using microarrays. D. pseu- 
doobscura and D. willistoni are particularly informative because 
they have independently derived X-autosome fusions (creating 
neo-X chromosome arms) that involved the same ancestral auto- 
some (see Meisel et al. 2009). We combined these expression 
measurements with available data on sex-biased and tissue-specific 
expression from both D. melanogaster and Mus musculus to more 
comprehensively examine how and why gene content differs be- 
tween the X chromosome and autosomes in Drosophila and 
mammals. 

Results 

X-Iinkage and sex-biased expression in Drosophila 

To confirm the results of microarray analyses (Parisi et al. 2003; 
Ranz et al. 2003; Sturgill et al. 2007), we measured sex-biased 

expression in whole flies using RNA-seq. As in the microarray ex- 
periments, our RNA-seq data show that the Drosophila X chro- 
mosome and D. pseudoobscura neo-X contain a dearth of genes 
with male-biased expression. We also find that genes with female- 
biased expression are over-represented on the X and neo-X chro- 
mosomes, and there is an excess of genes with male-biased ex- 
pression on the arm homologous to D. melanogaster chromosome 
2L (MuUer element B) in all species (Fig. lA), confirming patterns 
observed previously only in D. melanogaster (Parisi et al. 2003; Ranz 
et al. 2003). Most sex-biased expression is the result of genes 
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B. Head (mlcroarray or RNA-seq) 
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Figure 1. Frequency of genes with sex-biased expression on each chromosomal arm of Drosophila 
genomes. The percentages of genes on each chromosome arm that have female-biased (v\/hite), male- 
biased (black), or unbiased (gray) expression are indicated by the height of the bars. (Dashed lines) 
Genome-v\/ide percent of genes in each sex-bias class. (*) Significant deviations between observed and 
expected frequencies (as determined by a permutation test); (*) P < 0.05, (**) P < 0.01 , (***) P < 0.001 . 
Chromosome arms are sorted such that homologous arms are in the same order for each species vyith 
the dot chromosome (Muller element F) omitted: (Dmel) D. melanogaster; (Dpse) D. pseudoobscura; 
(Dwil) D. willistoni; (Dmoj) D. mojavensis. Expression was measured with RNA-seq using RNAfrom whole 
flies (Dmel, Dpse, and Dmoj) or abdomen (Dwil) (A), microarrays or RNA-seq from head (B), RNA-seq 
from thorax (C), and microarrays using RNA from third instar larva (D). In the case of the RNA-seq 
data, reads were aligned to the genome using BWA, and differential expression was determined using 
edgeR (see Methods). Species in panels A and B are sorted by increasing evolutionary distance from 
D. melanogaster. 



male-biased expression in head (Fig. IB; 
Supplemental Fig. 1). Interestingly, we 
detect an excess of genes with male-biased 
expression on the D. pseudoobscura neo-X 
based on our mlcroarray data (Fig. IB), 
but the RNA-seq data show the opposite 
pattern (Supplemental Fig. 1). In addi- 
tion, genes with male-biased expression 
in head are neither significantly over- nor 
under-represented on the D. willistoni X 
or neo-X (Fig. IB; Supplemental Fig. 1), 
and genes with male-biased expression 
in thorax may be under-represented on 
the D. willistoni X chromosome (Fig. IC; 
Supplemental Fig. 2). Therefore, unlike 
the pattern observed when reproductive 
tissues are included, there is not a con- 
sistent under-representation of X-linked 
genes with male-biased expression in non- 
reproductive tissues. On the other hand, 
genes with female-biased expression in 
head or thorax tend to be over-repre- 
sented on the X chromosome (Fig. 1B,C; 
Supplemental Figs. 1, 2), consistent with 
the pattern we observe when reproduc- 
tive tissues are included (Fig. lA). 

We next tested whether ontogenic 
changes in sex-specific selection pres- 
sures (Chippindale et al. 2001; Rice and 
Chippindale 2001; Gibson et al. 2002) 
affect the X-linkage of genes with sex- 
biased expression by measuring gene ex- 
pression in female and male wandering 
third instar larvae of D. melanogaster us- 
ing microarrays. We failed to detect a sig- 
nificant deficiency of X-linked genes with 



expressed in reproductive tissues (Parisi et al. 2003, 2004), so 
we measured expression in abdomen from male and female 
D. willistoni using RNA-seq to determine if the ancestral and neo-X 
chromosome arms have an excess/deficiency of genes with sex- 
biased expression. Genes with male-biased expression in D. willistoni 
abdomen are, indeed, under-represented, and genes with female- 
biased expression are over-represented on both arms of the X 
chromosome (Fig. lA). 

To test whether genes with sex-biased expression in non- 
reproductive tissues are over- or under-represented on the X 
chromosome, we measured expression in heads of males and fe- 
males of four Drosophila species using microarrays and RNA-seq. 
A previous analysis of RNA-seq data found an excess of X-linked 
genes with male-biased expression in D. melanogaster head (Chang 
et al. 2011), and we confirmed this pattern with our RNA-seq data 
(Supplemental Fig. 1). While we detect a similar trend in our 
mlcroarray data, the excess is not significant (Fig. IB). However, 
the X-linkage of genes with male-biased expression in non- 
reproductive tissues in other species is not consistent with the 
pattern in D. melanogaster. For example, genes with male-biased 
expression in head are not significantly over-represented on the 
D. mojavensis X chromosome when we analyze our mlcroarray data 
(Fig. IB), and they are under-represented on the X chromosome in 
our RNA-seq data (Supplemental Fig. 1). The ancestral X chromo- 
some in D. pseudoobscura is similarly depauperate in genes with 



male-biased expression in larvae (Fig. ID), 
suggesting that the factors responsible for the paucity of X-linked 
genes with male-biased expression do not on average act on larval- 
expressed genes. However, genes with male-biased expression 
in both larvae and adults may be under-represented on the X 
chromosome (Supplemental Fig. 3), indicating that the dearth of 
X-linked genes with male-biased expression may not be solely 
the result of genes specifically expressed in adults. We also 
find that genes with male-biased expression in larvae are over- 
represented on chromosome arm 2L, and genes with female- 
biased expression are over-represented on the X chromosome, 
consistent with the patterns observed in adults (Fig. 1). These 
trends are also observed when we consider genes that have the 
same sex-biased expression in both larvae and adults (Supple- 
mental Fig. 3). 

On the X-linl<age of Drosophila genes withi testis-biased 
expression 

We confirmed that D. melanogaster genes with enriched expression 
in testis relative to ovary are under-represented on the X chro- 
mosome (Supplemental Fig. 4; Sturgill et al. 2007; Vibranovski 
et al. 2009). Genes with testis-biased expression in mammals 
are identified by comparing across multiple tissues (e.g., Lercher 
et al. 2003). To parallel that analysis, we tested for a deficiency of 
X-linked D. melanogaster genes with testis-biased expression using 
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microarray measurements of expression from 14 adult tissues 
(ChlntapalU et al. 2007). 

First, we applied an approach in which the minimal ratio of 
testis expression level (microarray signal intensity) to the expres- 
sion level in other tissues is used as a measure of the degree of testis- 
biased expression (Mikhaylova and Nurminsky 2011), and we 
failed to detect a significant deficiency of X-linked genes with 
testis-biased expression (Supplemental Fig. 5). Next, we calculated 
T (Yanai et al. 2005), a metric of expression breadth that ranges 
from 0 (for broadly expressed genes) to 1 (for narrowly expressed 
genes). We applied t cutoffs to call genes as narrowly expressed in 
a single tissue or broadly expressed in multiple tissues, and we did 
not observe a significant deficiency of X-linked genes with testis- 
biased expression (Fig. 2), even at the highest t cutoffs (Supple- 
mental Fig. 6). This result is not because of low power to detect 
significant differences — genes with testis-biased expression are the 
largest class of narrowly expressed genes, which means that we 
have the greatest power to detect significant differences from the 
expectation. In addition, X-linked genes with testis-biased ex- 
pression have statistically indistinguishable t values when com- 
pared with autosomal genes with testis-biased expression (Fig. 3A; 
Supplemental Fig. 7), demonstrating that there is not a difference 
in the degree of testis-biased expression between X-linked and 
autosomal genes. 

Interestingly, genes with testis-biased expression that encode 
components of the sperm proteome (Dorus et al. 2006; Wasbrough 
et al. 2010) are under-represented on the X chromosome (Fig. 2; 
Supplemental Fig. 8). Therefore, while genes with enriched ex- 
pression in testis relative to ovary are under-represented on the X 
chromosome (Sturgill et al. 2007; Vibranovski et al. 2009), the only 
genes with testis-biased expression that are deficient on the X 
chromosome when additional tissues are considered are those that 
encode proteins incorporated in the ejaculate. 

Dosage compensation, sex-specific selection, and tlie X-Iinl<age 
of Drosophila genes with accessory-giand-biased expression 

We detect a significant deficiency of genes with accessory-gland- 
biased expression on the D. melanogaster X chromosome (Fig. 2; 
Supplemental Figs. 5, 6), consistent with the paucity of X-linked 
ACP genes (Wolfner et al. 1997; Swanson et al. 2001; Mueller et al. 
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Figure 2. X-linkage of D. me/onogoster broadly and narrowly expressed 
genes. (Barplots) Percent of D. melanogaster genes with narrow expression 
(t > 0.7) in each of four sex-limited tissues, with testis-biased expression 
and detectable in the sperm proteome (testis-SP), narrowly expressed in 
any of 14 tissues (narrow), narrowly expressed in one of 10 non-sex-lim- 
ited tissues (no sex) or broadly expressed genes (t < 0.4) that are X-linked. 
(Dashed lines) Percent of the entire genome that is X-linked. (*) Observed 
values that significantly differ from the expectation based on the size of the 
X chromosome as determined by a permutation test; (**) P < 0.01, (***) 
P< 0.001. 




X 2L 2R 3L 3R 
chromosome arm 

Figure 3. Distribution of t for each D. me/anogoster chromosome arm. 
(Boxes) Interquartile range; (horizontal line in the middle of the boxes) 
median value; (whiskers) 1.5x the interquartile range (outliers were 
omitted), (dashed line) genome-wide average. Significant differences 
between t for a given chromosome and the rest of the genome were 
assessed using a Mann-Whitney test; (*) P < 0.05; (**) P < 0.005. Distri- 
butions of T were calculated for (A) genes with testis-biased expression 
(t > 0.7), (6) genes with accessory-gland-biased expression (t > 0.7), and 
(C) all genes. 



2005). X-linked genes with accessory-gland-biased expression also 
tend to have broader expression profiles than autosomal genes 
with accessory-gland-biased expression (Fig. 3B; Supplemental 
Fig. 7). Therefore, genes with accessory-gland-biased expression 
are under-represented on the X chromosome, and those that are 
X-linked are not as narrowly expressed in accessory gland as are 
autosomal genes. This suggests that the exclusion of accessory- 
gland-expressed genes from the X chromosome may be responsi- 
ble for the dearth of X-linked genes with male-biased expression. 
However, if we exclude genes with accessory-gland-biased expres- 
sion, the X chromosome is still depauperate in genes with male- 
biased expression (Supplemental Fig. 9), demonstrating that the 
deficiency of X-linked accessory-gland-biased genes is not solely 
responsible. 

X-linked genes with accessory-gland-biased expression have 
lower expression levels in accessory gland than similarly narrowly 
expressed autosomal genes (Fig. 4A; Supplemental Fig. 10), sug- 
gesting that limits imposed by dosage compensation may be re- 
sponsible for the dearth of genes with accessory-gland-biased 
expression on the X chromosome (Vicoso and Charlesworth 
2009; Bachtrog et al. 2010). In comparison, the accessory-gland- 
expression levels of broadly expressed genes do not differ signifi- 
cantly between the X chromosome and the autosomes (Fig. 4A). 
If dosage limits in accessory gland were responsible for the paucity 
of X-linked genes with male-biased expression, we would expect 
to see the greatest deficiency of X-linked genes with male-biased 
expression among those genes that are highly expressed in acces- 
sory gland. Using sex-biased expression calls from our RNA-seq 
data, we found that genes with male-biased expression that are 
expressed in accessory gland at low, moderate, and high levels 
are all under-represented on the X chromosome (Fig. 4B). This 
suggests that dosage limits in accessory gland are not solely 
responsible for the deficiency of X-linked genes with male-biased 
expression. However, genes with male-biased expression that are 
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Figure 4. Accessory-gland expression level and X-linkage. {A) Boxplots 
show the distribution of the accessory gland expression level (S) for genes 
on each chromosome arm in D. melanogaster (chromosome 4 was omit- 
ted). (Boxes) Interquartile range; (horizontal line in the middle of the 
boxes) median value; (whiskers) 1.5x the interquartile range (outliers 
were omitted), (dashed line) genome-wide average. Genes were divided 
into those with accessory-gland-biased expression (t > 0.7) or those that 
are broadly expressed (t < 0.4). Significant differences in the expression 
level for a given chromosome and the rest of the genome were assessed 
using a Mann-Whitney test; (*) P< 0.05; (**) P< 0.005. (fi) D. melanogaster 
genes were divided into those with low expression in accessory gland 
(5 < 1 00), moderate expression in accessory gland (1 00 < S < 1 000), and 
highly expressed in accessory gland (S 2: 1000). Within each panel is 
shown the percent of genes with female-biased (F), male-biased (M), and 
unbiased (U) expression that are X-linked. Genes were assigned to sex- 
biased expression classes using RNA-seq data (BWA alignments and edgeR 
to call differential expression). (*) Observed values that significantly differ 
from the expectation based on the size of the X chromosome as de- 
termined by a permutation test; (*) P < 0.05, (**) P < 0.01 , (***) P < 0.001 . 

moderately or highly expressed in accessory gland are more under- 
represented on the X chromosome than genes with male-biased 
expression that are lowly expressed in accessory gland (P < 10"^, 
Fisher's exact test, FET). These results are robust to using different 
data to make the sex-biased expression calls (Supplemental Fig. 
11). We therefore conclude that dosage limits in the accessory 
gland are partially responsible for the dearth of X-linked genes 
with male-biased expression. 



Expression breadtli, sex-biased expression, and X-linl<age 
in Drosophila 

We confirmed that the D. melanogaster X chromosome has a dearth 
of narrowly expressed genes and an excess of broadly expressed 
genes (Fig. 2; Supplemental Figs. 5, 6), and we also found that 
X-linked genes tend to have broader expression profiles than 
autosomal genes (Fig. 3C). Additionally, the deficiency of narrowly 
expressed genes on the X chromosome remains when genes that 
are narrowly expressed in sex-limited tissues are excluded (Fig. 2; 
Supplemental Figs. 5, 6), indicating that genes with accessory- 
gland-biased expression are not solely responsible for the de- 
ficiency of X-linked narrowly expressed genes. 



If the paucity of X-linked genes with male-biased expression is 
a by-product of the correlation between expression breadth and 
sex bias (Meisel 2011; Mikhaylova and Nurminsky 2011), we do 
not expect to observe a deficiency of genes with male-biased ex- 
pression on the X chromosome when we control for expression 
breadth. However, there is a significant deficiency of X-linked 
genes in D. melanogaster with male-biased expression among both 
broadly and narrowly expressed genes (Fig. 5A; Supplemental Figs. 
12, 13), suggesting that the dearth of X-linked genes with male- 
biased expression is independent of expression breadth. The excess 
of X-linked genes with female-biased expression is also inde- 
pendent of expression breadth (Fig. 5A; Supplemental Figs. 12, 13), 
further demonstrating that the non-random X-linkage of genes 
with sex-biased expression is not a by-product of the correlation 
with expression breadth. 

On the other hand, narrowly expressed genes with non- 
sex-biased (unbiased) expression are not significantly under- 
represented on the X chromosome (Fig. 5; Supplemental Figs. 13, 
14). Additionally, there is not a significant difference in the pro- 
portion of narrowly expressed genes with unbiased expression that 
are X-linked when compared with the proportion of more broadly 
expressed genes with unbiased expression (P= 0. 102, FET) (Fig. SB), 
and this result is robust to most t cutoffs (Supplemental Table 1). 
We therefore conclude that the deficiency of X-linked narrowly 
expressed genes is an artifact of the correlation between expression 
breadth and sex bias, but the deficiency of X-linked genes with 
male-biased expression and excess with female-biased expression 
are independent of expression breadth. 

Reconciling tlie differences in gene content between 
tlie Drosophila and mouse X cliromosomes 

To examine the apparent differences in gene content between 
the Drosophila and mammalian X chromosomes (Vicoso and 
Charlesworth 2006; Gurbich and Bachtrog 2008), we analyzed ex- 
pression measurements from mouse in a similar manner as the 
Drosophila data described above. First, using microarray data (Yang 
et al. 2006), we confirmed that the mouse X chromosome has an 
excess of genes with female-biased expression in brain, liver, and 
muscle (Fig. 6A), much like the excess of X-linked genes with 
female-biased expression in Drosophila (Fig. 1). In addition, there 
is an excess of X-linked genes with male-biased expression in brain 
and a deficiency of X-linked genes with male-biased expression in 
muscle (Fig. 6A), demonstrating that in both Drosophila (Fig. 1) and 
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Figure 5. X-linkage, expression breadth, and sex-biased expression of 
D. melanogaster genes. (Barplots) Percent of D. melanogaster genes in 
various expression classes that are X-linked. (Dashed lines) Percent of all 
genes within that panel that are X-linked. (*) Observed values that sig- 
nificantly differ from the expectation based on the size of the X chromo- 
some as determined by a permutation test; (*) P< 0.05, (***) P< 0.001 . (A) 
Percent of broadly or narrowly expressed genes with either female-biased 
(F), male-biased (M), or unbiased (U) expression is plotted. (B) Percent of 
genes with unbiased expression that are either broadly (t < 0.4) or nar- 
rowly (t > 0.7) expressed is plotted. Genes were assigned to sex-biased 
expression classes using RNA-seq data (BWA alignments and edgeR to call 
differential expression). 



Genome Research 1259 



www.genome.org 



Meisel et al. 



6 

4 

T, 2 
$ 

X 

■5^ 6 
4 
2 
0 



adipose 

□■□ 




brain 


liver 


muscle 


■n 


U 





B 




3 13 cr 



CL CL CO ai 



U F 

sex-bias 



broad (t < 4) not narrow (t < 7) narrow (x > 7) 



cn cn ^ 
_^ 4i O- 



cr 

















In 




Id 




1 





M U 
sex-bias 



narrow broad 



Figure 6. Sex-biased expression, expression breadth, and X-linl<age in 
mouse. (Barplots) Percent of mouse genes in various expression classes 
that are X-linl<ed. (Dashed lines) Percent of the entire genome that is 
X-linked. (*) Observed values that significantly differfrom the expectation 
based on the size of the X chromosome as determined by a permutation 
test; (*) P < 0.05, (**) P < 0.01 , (***) P < 0.001 . (A) Percent of genes with 
female-biased (F), male-biased (M), and unbiased (U) expression that are 
X-linked is shown for data from four mouse tissues. (B) Percent of genes 
with narrow expression (t > 0.7) in each of four sex-limited tissues (white 
and black), with narrow expression in one of 25 tissues (gray), narrow 
expression in one of 21 non-sex-limited tissues (gray, no sex) or broad 
expression (t s 0.4; white) that are X-linked is plotted. Genes with testis- 
biased expression were identified using all testis EST libraries (testis), EST 
libraries from adult testis (UniGene IDs 251 1 and 2547), or EST libraries 
from 1 3-d and 1 5-d embryonic testis (UniGene IDs 2555 and 2594). (C) 
Percent of X-linked broadly expressed (t 0.4), not narrowly expressed (t 
:~0.7), or narrowly expressed (t>0.7) genes with either female-biased (F), 
male-biased (M), or unbiased (U) expression that are X-linked is plot- 
ted. (D) Percent of X-linked genes with unbiased expression that are ei- 
ther broadly or narrowly expressed is plotted. 



mammals, the X-linkage of genes with male-biased expression is 
tissue-dependent. 

We used mapped expressed sequence tags (ESTs) from adult 
mouse testes and unfertilized mouse ovaries to identify genes with 
enriched expression in testis relative to ovary. As in D. melanogaster 
(Supplemental Fig. 4; Sturgill et al. 2007; Vibranovski et al. 2009), 
there is a deficiency of X-linked genes with testis-enriched ex- 
pression (Supplemental Fig. 15). We next used mapped ESTs from 
25 mouse tissues to calculate t, and we found that there is an excess 
of X-linked genes with testis-biased expression (Fig. 6B; Supple- 
mental Fig. 16). However, the testis samples include a mixture of 
cell types from different developmental stages, including embryos 
(which should contain only pre-meiotic cells unaffected by MSCI) 
and adults (with pre-meiotic, meiotic, and post-meiotic cells) (Khil 
et al. 2004). To adequately mirror our analysis of the Drosophila 
data (Fig. 2), we only included EST libraries from adult testis, and 
we failed to detect a significant excess or deficiency of X-linked 
genes with testis-biased expression (Fig. 6B). This is true for most 
EST libraries and t cutoffs (Supplemental Figs. 16, 17). On the 
other hand, genes with biased expression in embryonic testis are 
over-represented on the X chromosome (Fig. 6B; Supplemental 
Fig. 16). Therefore, when similar analyses are performed on com- 
parable tissue samples, genes with testis-biased expression in adults 
are found at their expected frequencies on both the Drosophila 
and mouse X chromosomes. 



The mammalian prostate and Drosophila accessory gland are 
analogous male-specific organs that secrete seminal proteins nec- 
essary for male fertility. The D. melanogaster X chromosome has 
a paucity of genes with accessory-gland-biased expression (Fig. 2), 
and, similar to the pattern observed in humans (Lercher et al. 
2003), there is an excess of X-linked mouse genes with prostate- 
biased expression (Fig. 6B; Supplemental Fig. 16). We hypothesize 
as to why the patterns differ between Drosophila and mammals in 
the Discussion. 

While the Drosophila X chromosome appears to be deficient 
in narrowly expressed genes (Fig. 1; Mikhaylova and Nurminsky 
2011), the mammalian X chromosome is enriched for narrowly 
expressed genes (Fig. 6B; Supplemental Fig. 16; Lercher et al. 2003; 
Deng et al. 2011). Sex bias and expression breadth are correlated in 
mouse (Mank et al. 2008; Meisel 2011), so we tested whether the 
non-random X-linkage of genes with sex-biased or narrow ex- 
pression is an artifact of this correlation. When we exclude genes 
that are narrowly expressed in sex-limited tissues, there is no 
longer significant excess of X-linked narrowly expressed mouse 
genes (Fig. 6B; Supplemental Fig. 16). Additionally, narrowly 
expressed genes with non-sex-biased expression are not signifi- 
cantly over-represented on the mouse X chromosome (Fig. 6C,D; 
Supplemental Figs. 18, 19). Therefore, as in Drosophila, narrowly 
expressed mouse genes merely appear to be non-randomly X-Unked 
because of the correlation between sex bias and expression breadth. 
Lastly, the mouse X chromosome has an excess of genes with 
female-biased expression that are not narrowly expressed (Fig. 6C; 
Supplemental Fig. 18), much like the excess of X-linked Drosophila 
genes with female-biased expression (Fig. 5A). 

Discussion 

Our analyses combined measurements of sex-biased expression 
and expression breadth to determine more precisely what types 
of genes are over- or under-represented on the Drosophila and 
mouse X chromosomes so that we can infer the historical effects 
of sex-specific selection pressures. We have attempted to consider 
all possible confounding factors and correlations exhaustively, 
and our results suggest the following conclusions (summarized in 
Table 1). 

X-linkage and sex-biased expression in non-reproductive tissues 

Consistent with previous results (e.g., Parisi et al. 2003; Yang et al. 
2006; Chang et al. 2011), we showed that genes with sex-biased 
expression in non-reproductive tissues in both Drosophila and 
mammals are often non-randomly X-linked. Genes with female- 
biased expression are consistently over-represented on the X 
chromosome in both taxa (Figs. 1, 5A, 6), possibly because the 
female-biased transmission of the X chromosome favors the re- 
tention of alleles/genes with female-beneficial functions (Rice 
1984; Connallon and Clark 2010). The X-linkage of genes with 
male-biased expression varies across tissues, changes throughout 
development, differs between taxa, and can depend on the 
methodology used to measure expression (Figs. 1, 5A, 6; Supple- 
mental Fig. 1). The proportion of X-linked genes with male-biased 
expression in Drosophila head, for example, seems particularly 
labile (Fig. IB; Supplemental Fig. 1). We do, however, find some 
support for the hypothesis that stronger net selection in females 
prevents the accumulation of X-Unked genes with male-biased 
expression in non-reproductive tissues (Figs. 1, 5A; Supplemental 
Figs. 2, 9; Parisi et al. 2003). 
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Table 1. Factors responsible for non-random X-linkage 



Drosophila Mouse 



Expression Pattern Driven by Pattern Driven by 



Male-biased Tissue-dependent 

Female-biased X> A 

Narrow X<A 
Testis X = A 

Accessory gland or prostate X<A 



Sexual conflict in soma 

Acp genes on autosomes 

Sexual conflict in soma 

Ovary-biased on X 

Correlation with sex bias 

X < A from comparing testis and ovary 

Dosage limits 



Tissue-dependent Sexual conflict in soma 

Prostate-biased on X 
X> A Sexual conflict in soma 

X> A Correlation with sex bias 

X = A X > /lif exclude MSCI 

X> A Male-specific selection 



Interestingly, we detect an excess of X-linked genes with fe- 
male-biased expression in larvae (Fig. IC). If this excess is the result 
of female-beneficial substitutions being favored because of the 
X-chromosome's transmission bias, it would contradict the pre- 
diction that sexually antagonistic selection pressures are minimal 
in juveniles (Chippindale et al. 2001; Rice and Chippindale 2001; 
Gibson et al. 2002). However, third instar larvae have developing 
gonads (Hartenstein 1993), and this pattern could be driven by the 
excess of X-llnked genes with ovary-biased expression (Fig. 2). In 
addition, while genes with male-biased expression in larvae are 
found on the X chromosome at the expected proportion (Fig. IC), 
genes with male-biased expression in both larvae and adults may 
be under-represented on the X chromosome (Supplemental 
Fig. 3). The male-biased expression profile found in adults be- 
gins to manifest in the larval stages (Graveley et al. 2011), sug- 
gesting that selection on genes expressed in adults could affect 
the X-linkage of some larval-expressed genes. Additional data 
should be collected from larvae lacking gonadal tissues and from 
earlier stages of development to further interrogate how ontogenic 
changes in sexually antagonistic selection can affect the X-linkage 
of genes encoding proteins that perform sex-beneficial functions. 

X-iini(age and biased expression in male reproductive tissues 

The D. melanogaster X chromosome has a deficiency of ACP genes 
(Wolfner et al. 1997; Swanson et al. 2001; Mueller et al. 2005) and 
genes with accessory-gland-biased expression (Fig. 2; Mikhaylova 
and Nurminsky 2011). We hypothesize that this is the result of 
dosage limits on X-linked genes in males, and these dosage limits 
in accessory gland appear to be partially responsibly for the paucity 
of genes with male-biased expression on the Drosophila X chro- 
mosome (Fig. 4B). Our expression data from larvae support this 
hypothesis. Male third instar larvae contain gonads with cells that 
have begun spermatogenesis (Fuller 1993), but they have yet to 
develop accessory glands (Hartenstein 1993). There is not a signifi- 
cant deficiency of X-linked genes with male-biased expression in 
larvae (Fig. ID), suggesting that genes expressed in accessory gland 
make an Important contribution to the dearth of X-Iinked genes 
with male-biased expression in adults. Genes with testls-biased 
expression that encode components of the sperm proteome are 
also under-represented on the Drosophila X chromosome (Fig. 2). 
Both ACPs and sperm proteins must be synthesized at high levels 
(Wolfner et al. 2005), which could exacerbate the dosage limits 
imposed by being X-Iinked. 

The Drosophila accessory gland plays an important role in 
male fertility, secreting proteins that are incorporated into the 
seminal fluid, which mediate post-mating inter- and intrasexual 
conflict (Swanson and Vacquier 2002; Ravi Ram and Wolfner 
2007). While it is tempting to Invoke sexually antagonistic selec- 



tion models (e.g.. Rice 1984; Patten and Haig 2009; Connallon and 
Clark 2010) to explain the paucity of X-linked ACP genes, these 
models assume intralocus sexual antagonism. Drosophila male re- 
productive proteins, such as ACPs, are thought to mediate inter- 
sexual conflict by interacting with proteins expressed in the 
female reproductive tract (e.g., Yapici et al. 2008). Therefore, new 
models will be required to examine whether interlocus sexual 
conflict can explain the non-random X-linkage of male repro- 
ductive proteins. 

Mammalian genes with prostate-biased expression, unlike 
Drosophila ACPs, are over-represented on the X chromosome (Fig. 6B; 
Letcher et al. 2003), as are mouse genes with testis-biased expression 
prior to meiosis (Fig. 6B; Khil et al. 2004). X-chromosome dosage in 
eutherian mammals is equilibrated by up-regulation of X-linked ex- 
pression m both sexes followed by random mactivation of one X in 
each female cell (Payer and Lee 2008), which does not appear to 
impose dosage limits on highly expressed X-linked genes in males. 

The excess of genes with prostate-biased and early sper- 
matogenic expression on the mammalian X chromosome (Fig. 6B; 
Lercher et al. 2003; Khil et al. 2004) could be the result of sex- 
specific selection pressures. Unlike Drosophila ACPs, mammalian 
seminal proteins do not alter female physiology at a large scale 
(Pitnlck et sil. 2009), but they do mediate male-male sexual selec- 
tion (Pizzari and Parker 2009). Therefore, the selection pressures 
on genes with prostate- and testis-biased expression likely act in a 
male-limited manner. If this is the case, the female-biased trans- 
mission of the X chromosome will not affect the X-linkage of these 
genes. We hypothesize that the hemlzygous transmission of 
the X chromosome in males allows for the selective fixation of 
male-beneficial recessive mutations in X-linked mammalian genes 
(Charlesworth et al. 1987), promoting the accumulation of genes 
expressed in male reproductive tissues on the X chromosome. 

MSCI and tlie X-linl<age of genes witli testis-biased expression 

We detect a paucity of X-Iinked D. melanogaster genes with en- 
riched expression in testis relative to ovary (Supplemental Fig. 4), 
but not when we examine genes with testis-biased expression 
compared with multiple tissues (Fig. 2). The previously noted de- 
ficiency of X-linked genes with testis-enriched expression (Sturgill 
et al. 2007; Vibranovski et al. 2009) is likely a methodological 
artifact for the following reasons. First, genes with ovary-biased 
expression are over-represented on the X chromosome (Fig. 2). 
Second, testis-expressed genes tend to be narrowly expressed (Meisel 
2011), and narrowly expressed genes are under-represented on 
the X chromosome (Fig. 2). Furthermore, ovary-expressed genes are 
broadly expressed (Meisel 2011), and there is an excess of X-linked 
broadly expressed genes (Fig. 3C). While the relationship between 
expression breadth and X-Iinkage appears to be an artifact of the 
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correlation between expression breadth and sex-biased expression 
(Fig. 5), direct comparisons between testis and ovary expression do 
not control for this. Therefore, by comparing testis and ovary ex- 
pression, one artificially obtains evidence for a paucity of X-linked 
testis-enriched genes. 

The lack of a deficiency of X-linked genes with testls-biased 
expression does not, however, reveal Insights Into the debate over 
MSCI In Drosophila (Hense et al. 2007; Vlbranovskl et al. 2009; 
Melklejohn et al. 201 1; Mikhaylova and Nurmlnsky 201 1). Mouse 
genes with testls-biased expression are also not under-represented 
on the X chromosome even when meiotic tissues are included 
(Fig. 6B; Khil et al. 2004), and MSCI is well documented in mouse 
(Turner 2007). On the other hand, our results as well as others — 
e.g., the X chromosome Is deficient In genes with mcile-biased 
expression even when testis-expressed genes are excluded (Sturgill 
et al. 2007) — suggest that, even if MSCI occurs In Drosophila, it 
is not an important proximate cause of the paucity of X-llnked 
genes with male-biased expression. 

Conclusions 

In summary, our results suggest that most of the observed differ- 
ences In X-chromosome content between Drosophila and mam- 
mals are confounded by the correlation between sex bias and 
expression breadth, or they are affected by taxon-spedflc data 
collection and approaches to the analysis. The female-biased 
transmission of the X chromosome appears to favor the accumu- 
lation of female-beneficial alleles/genes on the X chromosome in 
both taxa, but the effects on male-beneficial alleles/genes are not 
consistent between taxa or across tissues. While models assuming 
Intralocus sexual antagonism may explain many of these patterns, 
further work is required to develop theoretical models of the effects 
of interlocus conflict on the X-linkage of genes that are under sex- 
speclflc selection pressures. 

Methods 

Samples 

Flies were grown under the following conditions: D. melanogaster 
(yi iv'''^ for whole fly, OreR for head and larvae) at 22°C on corn- 
meal media; D. pseudoobscura (Drosophila Species Stock Center 
[DSSC] 14011-0121.94) at 22°C on banana-opuntia media; 
D. wilHstoni (DSSC 14030-0811.25) at 18°C on cornmeal media; 
D. mojavensis (DSSC 15081-1352.22) at 22°C on banana-opuntia 
media. For the whole flies from D. melanogaster, D. pseudoobscura, 
and D. mojavensis, 5- to 7-d-old females cind males were sampled. 
For the head samples, newly eclosed D. melanogaster, D. pseu- 
doobscura, and D. mojavensis were aged in vials (males and females 
together) for 7 d. The sexes were separated on day 7 under CO2 
and stored in separate vials for 24 h. Flies were flash-frozen on dry 
Ice after 24 h and stored at — 80°C; heads were dissected from 
frozen flies using forceps on ceramic plates chilled on dry Ice. 
Males and females of newly eclosed D. willistoni were separated 
under CO2 and aged for 5-7 d; whole flies were placed live In 
Ringers solution, and heads, thoraxes, and abdomens were dis- 
sected; samples from each segment were flash-frozen In liquid N2 
and stored at -80°C. D. melanogaster wandering third instar larva 
were collected and sexed based on the size of the visible gonad 
(Demerec 1950). Males and females were separately flash-frozen 
in liquid N2 and stored at -80°C. RNA was extracted from the 
samples using TRlzol (Invltrogen) following the manufacturer's 
Instructions. 



Microarrays 

Mlcroarray measurements of gene expression were performed on 
RNA samples from the heads of D. melanogaster using the stock 
Nlmblegen 60-mer array (Dmel r. 4.2.1), heads of D. pseudoobscura 
and D. mojavensis using custom Nlmblegen 50-mer arrays (Zhang 
et al. 2007), and larvae of D. melanogaster using Agilent 4-sector 60- 
mer arrays (G2519F). Four biological replicates were collected of 
the head mlcroarray data, and three biological replicates were 
collected of the larval mlcroarray data. Total RNA was reverse- 
transcribed to cDNA, samples were dye-labeled, and hybridization 
was performed per the manufacturers' instructions. For the head 
data, RMA normalization was performed using Nimblescan v 2.4 
software (Roche Nlmblegen Inc.). We tested the null hypothesis 
that a gene is expressed at the same level in females and males 
using a moderated t-test Implemented In the LIMMA package of 
Bloconductor In the R statistical software environment (R De- 
velopment Core Team 2009) with the empirical Bayes function to 
pool sample variances toward a common value (Smyth 2004). For 
the larval data, the LIMMA package was used to normalize channel 
intensities between arrays (Smyth and Speed 2003) and test for 
differential expression between samples (Smyth 2004). P- values 
from statistical tests were corrected for the false discovery rate 
(FDR) using the Benjamini and Hochberg (1995) adjustment. Ad- 
ditional information can be found in GEO accessions GSE23309 
and GSE3 1722. 



RNA-seq 

lUumlna libraries were prepared using RNA extracted from males 
and femcdes from the following samples: D. melanogaster whole fly, 
D. melanogaster head, D. pseudoobscura whole fly, D. pseudoobscura 
head, D. willistoni head, D. willistoni thorax, D. willistoni abdomen, 
D. mojavensis whole fly, D. mojavensis head. Standard lUumina li- 
brary sample preparation procedures were used: Poly-adenylated 
mRNA was Isolated using either the QIAGEN Oligotex mRNA kit or 
ollgo-dT Dynabeads (Invltrogen); mRNA was fragmented using 
alkaline hydrolysis; double-stranded cDNA was created using ran- 
dom hexamer primers; sequencing adaptors were llgated to the 
cDNA; fragments were size-selected by gel purification; and the 
size-selected fragments were amplified by PCR. Samples were run 
on the following instrumentation: lUumina Genome Analyzer 1 
(one lane each of male and female D. melanogaster, D. pseu- 
doobscura, and D. mojavensis head samples), lllumina Genome 
Analyzer 11 (one lane each of male and female D. melanogaster, 
D. pseudoobscura, and D. mojavensis whole fly samples; three 
lanes each of male and female D. melanogaster head samples; three 
lanes each of male and female D. pseudoobscura head samples; four 
lanes each of male and female D. mojavensis head samples; one lane 
each of male and female D. melanogaster, D. pseudoobscura, and 
D. mojavensis whole fly samples; one lane each of male and female 
D. willistoni head, thorax, and abdomen samples), lllumina HiSeq 
(one lane each of male and female D. melanogaster, D. pseu- 
doobscura, and D. mojavensis whole fly samples; one lane each of 
male and female D. willistoni head samples). Additional details can 
be found In GEO accessions GSE20348, GSE20882, GSE19989, 
GSE28078, and GSE31723 (Graveley et al. 2011; Malone and 
Oliver 2011). A full list of the NCBI Sequence Read Archive runs 
can be found In Supplemental Table 2. 

Reads were processed using the standard lllumina pipeline 
(see GEO accessions), and all reads longer than 36 bp were short- 
ened to 36 bp to avoid biases introduced by different read lengths 
and to eliminate low-quality sequences at the 3' end of longer 
reads. For paired-end runs, we only used the reads In the forward 
direction. The reads were aligned to the annotated proteln-codlng 
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sequences using BWA (Li and Durbin 2009) and Bowtie (Langmead 
et al. 2009). The following annotations were used: D. melanogaster 
r5.36, D. pseudoobscura r2.19, D. willistoni rl.3, and D. mojavensis 
rl.3. Only the longest annotated protein-coding sequence was 
used for each gene. Reads that mapped to multiple locations were 
not counted, and genes with fewer than 50 mapped reads (males 
and females combined) were excluded from subsequent analyses. 
We then tested for differential expression between male and fe- 
male samples using DESeq (Anders and Huber 2010), edgeR 
(Robinson et al. 2010), and Cufflinks (Trapnell et al. 2010), using 
P-values corrected for cin FDR of 0.05 (Benjamini and Hochberg 
1995). Our results are robust to the alignment algorithm and sta- 
tistical analysis used (Supplemental Fig. 1) unless otherwise men- 
tioned above. 

Additional sex-biased expression data sets 

We obtained microarray measurements of sex-biased expression in 
whole D. melanogaster and head from SEBIDA (Gnad and Parsch 
2006). We also obtained microarray measurements of sex-biased 
gene expression from four mouse tissue samples (Yang et al. 2006). 
For the mouse data, we used the results of a Mann- Whitney test to 
assign genes to sex-biased expression categories (Mank et al. 2008), 
and we applied an FDR cutoff of 0.05. This was done separately for 
each tissue. We called genes as sex-biased in mouse (Fig. 6C,D) if 
they were assigned to a sex-biased expression class in one tissue 
and were either put in the same sex-biased expression class or not 
called as sex-biased in all other tissues. 

Expression breadtFi 

Microarray signal intensities (S) from 14 adult D. melanogaster tis- 
sues (brain, eye, thoracicoabdominal ganglion, salivary gland, 
crop, midgut, tubule, hindgut, heart, fat body, ovary, testis, ac- 
cessory gland, and spermatheca) were obtained from FlyAtlas 
(Chintapalli et al. 2007). Measurements from spermatheca of 
virgin and mated females were averaged to obtain a single sper- 
matheca signal because they are well correlated (Meisel 2011). 
Genes were said to be moderately expressed in a tissue if S > 100 
and highly expressed if S a 1000, per the curators' instructions 
(Chintapalli et al. 2007). 

The number of ESTs mapping to each gene in the mouse 
genome was obtained for 25 tissues (adipose, bladder, blood, bone, 
bone marrow, brain, dorsal root ganglion, eye, vagina, intestine, 
heart, inner ear, liver, lung, lymph node, mammary gland, muscle, 
pancreas, prostate, skin, spinal cord, spleen, sympathetic ganglion, 
testis, and thymus) from UniGene build 190. Additionally, the 
following individual EST libraries were extracted from adult or 
embryonic testis: 254 (Yuan et al. 1995), 483, 1336 (Marra et al. 
1999), 1784, 2511, 2555, 2594, 18007 (Okazaki et al. 2002), 2547 
(Caminci et al. 2003), and 7009 (Strausberg et al. 2002). 

Expression breadth was calculated two ways. In the first ap- 
proach (applied only to the D. melanogaster data), we calculated the 
minimal ratio of the expression level in a particular tissue relative 
to all other tissues (Mikhaylova and Nurminsky 2011). In the 
second approach, we calculated a metric of expression breadth as 
follows: 

yN _ logS, 

_ lOg^nm 

^~ N-1 

where N is the number of tissues analyzed, S, is the expression level 
in tissue / (signal intensity in FlyAtlas and number of ESTs mapped 
to the gene standardized by the total number of mapped ESTs in 
that tissue for the UniGene data), and Smax is the maximum ex- 
pression level across all tissues (Yanai et al. 2005; Larracuente et al. 



2008; Mank et al. 2008; Meisel 201 1). In the FlyAtlas data, all S,- < 1 
were set to 1 so that log(Sj) a 0. If multiple probes were present for 
a given gene, the median S across all probes was used for that 
gene. In the UniGene data, all genes with fewer than three 

mapped ESTs across all 25 tissues were removed, and S; < 2 were 
set to 2 for the remaining genes. The t value for each gene was 
used to determine if the gene is narrowly or broadly expressed 
using various cutoffs. 

Sperm proteome 

Genes were said to encode proteins in the sperm proteome if they 
were detected in at least one of two mass spectrometry analyses of 
D. melanogaster sperm proteins (Dorus et al. 2006; Wasbrough et al. 
2010). 

Statistical analysis 

We permuted the chromosomal locations of the genes in our data 
sets 1000 times to estimate a null distribution of the number of 
genes in each expression class on each chromosome arm in Dro- 
sophila and each chromosome in mouse. Depending on the com- 
parison being made, this was either done for all genes or a subset of 
genes under consideration. All permutations were performed in 
the R statistical package (R Development Core Team 2009). 

Data access 

All microarray and RNA-seq data presented in this mcmuscript 
have been submitted to the NCBI Gene Expression Omnibus 
(GEO) (http://www.ncbi.nlm.nih.gov/geo/) under accession num- 
bers GSE23309, GSE31722, and GSE31723. 
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